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Abstract 



We consider a mass conserved Allen-Cahn equation ut = Au + 
e~'^{f{u) — eA(i)) in a bounded domain with no flux boundary condi- 
tion, where £\{t) is the average of f{u{-,t)) and — / is the derivative 
of a double equal well potential. Given a smooth hypersurface 70 con- 
tained in the domain, we show that the solution with appropriate 
initial data approaches, as e \ 0, to a limit which takes only two val- 
ues, with the jump occurring at the hypersurface obtained from the 
volume preserving mean curvature flow starting from 79. 

1 Introduction. 

In this paper, we study the hmit, as e — >■ 0, of the solution to the mass 
conserved Allen-Cahn equation (P^) 
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on on X 

on i7 X {0}, 




where 

/ /(^')= / f{u%x,t))dx. 
J n Jn 

Here f2 is a smooth bounded domain inffi.*^ {n>l), di, the outward normal 
derivative to dVL, and —f{u) is the derivative of a smooth double equal well 
potential; more precisely, 

feC^(K),fi±i) = o,fi±i)<o,iy = £f<o Vwe(-i,i). (2) 

A typical example is f{u) = u ~ u^. The initial data satisfies, for some 
smooth hypersurface 70 CC ^2, 

lim/(x) = |7j Vx6n\7o. (3) 

£-^0^ ^ ' [ +1 outside 7o ^ ^ 

Problem ([1]) was proposed, along with its well-posedness, by Rubinstein and 
Sternberg [19] as a model for phase separation in binary mixture. The model 
is mass preserving and energy decreasing since 

d r 

Vt > 0, — / u'(x,t)dx = 

and 

/u 
f{s)ds is the double equal well potential. 

Formally, one can show that, as £ — 0, the solution to ([T]) and tends 
to a limit 

where F := [J (74 x {t}) is the solution to the volume preserving mean cur- 
vature motion equation 

V={n- 1)K^, - ^^^r^ I K^.dH''-^ on 7^ (5) 
\lt\ ht 
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starting from 70. Here V is the normal velocity of '-ft (positive when 7^ is 
shrinking) and K^^ the mean curvature (positive at points where jt is locally 
the boundary of a convex domain). 

The local in time existence of a unique smooth solution to has been first 
established in a two-dimensional setting in [11]. The general result in arbi- 
trary space dimension is obtained in [13], where the large time behaviour of 
solutions for initial data close to a sphere was also investigated. When the 
initial data is convex, it is shown in [16] that ([5]) admits a unique global in 
time convex solution. Related properties of other volume-preserving curva- 
ture driven flows are established in [T2] . 

Concerning the connection between ([T]) and Bronsard and Stoth [3] con- 
sidered a radially symmetric case with multiple interfaces (rings) and proved 
dl]). Let us also mention [T^ where a similar result is established for a 
different nonlocal mass conserved AUen-Cahn equation, using the method 
introduced in p]. In the present paper, we shall consider general smooth 
initial interfaces 70 CC f2 and prove the following: 

Theorem 1 Let T = Uo<t<T(7t ^ {t}) be a smooth solution to ^ satisfying 
7t CC for all t G [0, T]. Then there exists a family of continuous functions 
{(7^}o<e<i such that the solution to (Qj) satisfies ^ for all t G [0,T]. 

For the Allen-Cahn equation ul = Am^ — 6^'^f{u'^), (jlj) holds with F being 
the solution to the motion by mean curvature flow V = {n — 1)K^^. A simple 
method to verify this is to use a comparison principle and construct sub- 
super solutions [H [H] . There are different notions of weak solutions such as 
viscosity [11] and varifold [17] which can be used to establish the global in 
time limit. Nevertheless, ([1]) does not have a comparison principle (due to the 
volume preserving property) and the simple method does not seem to work. 
Here we shall employ a method first used by de Mottoni and Schatzman [10] 
for the Allen-Cahn equation, and later on by Alikakos, Bates, Chen [1] for 
the Cahn-Hillard equation and Caginalp and Chen [6] for the phase field 
system. 

Namely we first rewrite the equation for in Problem (P^) as 

< = Au' + e-\f{u') - eX,{t)) in n x 1R+, (6) 
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where we define 



Vt>0, X^(t) = yjiu'{.,t)). (7) 
The basic strategy of the proof goes as follows. 



1. For a large enough A; G N, construct an approximate solution A|) 
satisfying 

ul, - Aul - e-\f{ul) - e\%) = 5% in := ^ [0, T], 
^ ul 4x = Vt G [0, T], d^ul = on := dVt x [0, T] 



where 51 = 0{l)e^. Note that, by integration, 

J n 



2. For each t G [0,T] and small positive e, estimate the lower bound of 
the spectrum of the self-adjoint operator —A — e~^f'{u%{-, t)); namely, 
show that for some positive constant C*, 

inf inf inf f - e-'f'{ul{.,t))<p') > -C*. (9) 

3. Set R = u'^ — u% and show that R tends to as £ ^ 0. 

The organization of this paper is as follows. In section 2, we present an error 
estimate required in step 3. In section 3, we recall a known spectrum estimate 
[9l [5] that can be adapted here to prove step 2 in the strategy described 
above. After some geometrical preliminary computations in section 4, we 
finally construct the approximate solutions in section 5. 

2 Error Estimate 

The error estimate relies on the following result which is proved in the ap- 
pendix. 
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Lemma 1 Lei f2 C IR." with n > 1) be a bounded domain, letp = min{4 /n, 1 } . 
Then there exists C = Cn(^) > such that for every R G H^{Q) with 
J^Rdx = 0, 

\\RfX<C\\R\\l4VR\\l., (10) 
where for any q > 1, L'^ = L''{Q). 

Rubinstein- Sternberg established in [19] L°° bounds for the solution u*^ to 
Problem (P^) using invariant rectangles. Therefore we can modify / outside 
of a compact interval and assume for simplicity that 

lim f(u) = =Foo 

M— >±00 

and that there exists M > such that 

V|m| > M, uf"{u) < 0. 

Since p G (0, 1], for any Cq > 0, there exists C = C{Co,p) such that for all 
\u\ < Co and R eR, 

{f{u + R)- f{u) - f'iu)R)R < C\Rr'. 

Indeed, note that for i? in a compact interval, 

{fiu +R)- f{u) - nu)R)R = f:>±MRS < c\Rr', 

whereas for \R\ — > +oo, f(u + R)R — oo, uniformly in \u\ < Cq so that 
{f{u + R)- f{u) - f\u)R)R < {-f{u) - f\u)R)R < CR^ < C\R\p+^ . 

Lemma 2 Assume that k > max{4,n} and {'u|}o<e<i satisfies ^ and ^ 
with 

Let {m^}o<£<i be solutions to (QP with initial data {g"^} satisfying 

/(•) = «U-,O) + 0^(-), /n0^ = O, U'\\LHn)<e'. 
Then for all sufficiently small positive e, 

sup \\u'{-,t)-ul{-ML'm<C{T)e\ 

0<t<T 
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Remark 1 By a bootstrap argument, one can show that other norms of {u'^ — 
u%) tend to as £ \ 0. 

Proof. In the sequel, C denotes a generic positive constant independent of e. 
Set p = min{l, 4/n} and R = u"^ —u%. Then f^ R{x, t)dx = for all t e [0, T]. 
Also, 

R{f{u^) - f{ul) - f'{ul)R} < C\R\'^^. 

Multiplying by R the difference of the equations for and m| and integrating 
the resulting equation over Q gives, after integration by parts. 



1 d 
2dt^ 



\R\\h+ f {\VR\'-e-'f{u^^)\R\'}< f {Ce-'\R\'+^+\R5l\}. 
By©, 

/ {|Vi?r - e-'f'iul)R'} = [ +(1 - e') [ > e^VR\\h - C\\R\\h. 
The interpolation (ITOl) then yields 

\j^\\R\\h < C\\61\\l4R\\l2 + C\\R\\l. - II Vi?||i.{^' - C,e-'\\R\\l,}. (11) 
We define 

:= sup{t G [0,T] | \\R{-,t)\\v2 < e^l^'C^^'" for all r G [0,t]}. 
Since /c > max{4, n\ = A/p, it follows that 

\\R{;0)\\L.<e' <e'/^C^'^' 
for e > small enough. Therefore, > 0. Also, from fill I) , we have for all 

te{o,T% 

±\\R\\,,<cm\v^ + \m\L^) 

The Gronwall's inequality then provides 

sup m-ML^ < e^^[||i?(-,0)|U. + rmh^dt] < C{T)e^ < \e"^C^^'^ 
if e is small enough. Thus, we must have = T. This completes the proof. 
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3 The linearized operator 
3.1 A Spectrum Estimate 

Assume that / satisfies ([2]). Then there is a unique solution 9q{-) :]R. ^ (0, 1) 
to 

e'^ + f{eo)=OonR, ^o(±oo) = ±l, 0o(O) = O. (12) 



The solution satisfies, for a = mm{y —f'{l), y — /'(— 1)}, 

^ 1) = ©(e-"!"!) as ± p oo, Vm gN. 
Let Oi eC^(K)n be any function satisfying 



'2 f"^%)e^ =0. (13) 



Let fi^ CC be a subset with boundary 7 = dfl^ . Denote by d{x) the 
signed distance (negative in Q~) from a; to 7 and by s{x), for x close to 7, 
the projection from x on 7 along the normal to 7. 

We look for the spectrum of the linearized operator of —Au—e^'^f'{u) around 
u = ip'^ given by 

, J eo{^)+ef{s{x))e,{^) + 0{l)e' if \d{x)\ < ^S, 
'?/' (x) = < 6 e (14) 

i ±1 + 0(l)e if ±d{x) > ^. 

The following spectrum estimate was first proven by de Mottoni and Schatz- 
man |9], then by Chen [5j in a more general situation that can be used in 

[DE]. 

Proposition 1 Let 7 G C^, and ]f and 0(1) in ( (771 ) hounded indepen- 
dently of 6. Then there exists a positive constant C* depending on ||7||c3; 
IIp^IIloo and ||0(1)||loo such that for every e G (0,1] and (p G H^iVL), 



7 



We define the linearized operator around 6'o(p) acting on v = v{p) by 

Cv:= -v" - f'ieo)v (15) 

In our application, 6i is the unique solution to 

C9i = l-a9'Q in R 

6,(0) = 0, a:=2/^6',\ 

Integrating Oq C 9i over ffi, and by parts, one can verify that (fT3ll is satisfied; 
see Piail]. 

We remark that the distance function d in f|T^ can be replaced by a "quasi- 
distance" function d'^ given by 

d^{x) = d{x) - ehi{s{x)) - e'^h2{s{x)) + 0{l)e^ 

where hi and h2 are smooth functions on s G 7. 

3.2 Solvability Condition 

Lemma 3 Assume that f satisfies Let 6q be the solution to [T^) . a = 



min{y— /'(I), y— /'(— 1)} and C be defined in 173]) . Assume that a function 
h{p,s,t) satisfies, as p ±00, 

D^D:D\[h{p,s,t) - h^it)] = 0(|pre-"l''l) 

for some i > and all (m, n, I) G and (s, t) in U x [0, T] . Then 

CQ = h{-,s,t) in R, Q(0,s,t) = 
has a unique bounded solution Q{p,s,t) if and only if 

V(s,t) G f/ X [0,T], / h{p,s,t)e'o{p)dp = 0. 

// the solution exists, then it satisfies, for all (m, n, /) G and {s,t) G 
Ux[0,T], 

D^D:Dl[Q{p,s,t) + = Odpre-I^l) as p ^ ±00. 

Proof Since C6q = 0, the ode CQ = h can be solved explicitly. We omit 
the details of the proof; see [9]. 
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4 Differential Geometry: local coordinates 

4.1 Parametrization around the limit interface 

Let r = Utg[o,T] It X {t} C Qt be the smooth solution to ([5]) on [0,T] and 
Q^{t) the two domains separated by 7^, with 7^ = 9fi~(t). For each fixed t, 
we use d(x, t) to denote the signed distance from x to 'jt (positive in Q^(t)). 
Then d(-, ■) is smooth in a tubular neighborhood of the interface. Locally we 
choose a parametrization of 7t by Xo(s,t) with s eU Cffi."""^ so that 

-) (16) 



9s 1 ds. 



n-l 



is a basis of the tangent space to 7t at Xo(s, t), for each s E U. We denote 
by n(s,t) the unit normal vector to 74, pointing towards Q^{t) so that 

n(s,t) = Vd(Xo(s,t),t). 

Up to a suitable multiplication factor si — Asi, we may assume that 

Next for each fixed t G [0,T], a local parametrization by coordinates (s,r) G 
[/ X (—35, 35) is obtained by 

X = Xo{s,t) +rn{s,t) = X{r,s,t), (18) 

which defines a local diffeomorphism from (—35, 35) x U onto the tubular 
neighborhood of 7^, 

F3*5 = {a;Gfi, |d(a;,t)| <35}. (19) 
We denote the inverse by 

r = d(x, t), s = S(x, t) = (5^(x, t), S\x, t), ^"-^(x, t)). (20) 
In particular, since for all fixed s G f/, t G [0, T] and for all r G (—35, 35), 

d{Xo{s, t) +r n(s, t), t) = r. 
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it follows by differentiation with respect to r that for all r G (—35, 3(5), 
Vd{Xo{s, t)+r n(s, t), t).n(s, t) = 1. 

Using that 

\Vd{x,t) \ = 1 for X close to 7^, (21) 
this equality imposes that for all (r, s) G {—36, 36) x U, 

Vd{Xo{s, t)+r n(s, t),t) = n(s, t) (22) 

proving that Vd is constant along the normal lines to 74. Thus the projection 
from X on 7^ is defined by 

Xo(S(x, t), t) = X - d(x, t) Vd(x, t). (23) 

It follows also from that for alH = 1, n and for x G Vg*^, 

" d'^d dd 

Thus the symmetric matrix D'^d{x, t) has eigenvalues ■ ■ ■ , 0} with 
unit eigenvectors {ri, ■ ■ ■ , r„_i, Vc/} forming an orthonormal basis of ffi." for 
X G ^3*^. In particular, for x & 'jt, the are the principal directions and the 
Ki are the principal curvatures of 7^. Note that {ri, ■ ■ ■ , r„_i} form a basis 
of the tangent hyperplane to 7^ at a; = Xo(s,t). By definition, K and /T^^ 
are respectively the sum of principal curvatures and the mean curvature of 
7t, given by 

n-l 

K = {n- l)K^, = Ad{Xois, t), t) = E «^^(^' (25) 

i=l 



Note that using flMl) . for x G 7(, we have that 

Vd-VAd = vA(^^)_vf^ 

^ ^9xj dxidxj ' ^ ^dxidxj 
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We denote 

n-l 

bis, t) = - Vd ■ VAd|xo{s,i),i = E (26) 

i=l 

Let V{s,t) be the normal velocity of the interface at the point Xo(s,t) so 
that using ([22]), 

V{s,t) = Xot{s,t).n{s,t) 

= Xot{s,t).Vd{Xo{s,t)+rn{s,t),t) = -dt{X{r,s,t),t) (27) 

where the last equality follows from differentiating with respect to t the 
identity 

d{Xois,t) +rn{s,t),t) = r. 

It follows that dt{x,t) is independent of r = d{x,t) for \r\ small enough. 
Changing coordinates from {x,t) to (r, s,t), we associate to any function 
t) the function 

4>{r, s, t) = </.(Xo(s, t) + rn(s, t), t) (28) 

or equivalently 

0(x, t) = (f>{d{x, t), S{x, t),t). 
By differentiation we obtain the following formulas 

dt(j) = {-Vdr + dl)4> 

A0 = {drr + Addr + A^)4> (29) 

with 

n-l 
i=l 

n-l 

= (Ev5%)0 

j=i 

n— 1 n— 1 

A^^ = (EA5%+ E V^'-V^'^.M.)0 (30) 

i=l iii=l 
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where VS\ Sj, Ad, d t are evaluated at x = X{r,s,t) and are viewed as 
functions of {r,s,t). Note that the mixed derivatives of the form d'^^j(f) do 
not appear eventually in fl29|) because for all j = 1, 2, n — 1, 

VS^{x,t).Vd{x,t) = 

(This follows from differentiating with respect to r the identity 

Vr G (-3(5,35), S^{Xo{s,t) +rn{s,t),t) = 

which holds for all fixed s G f/, t G [0, T] and j = 1, 2, n — 1.) 



4.2 The stretched variable 

Following the method used in [7j, we now define the stretched variable p by 
considering a graph over 7^ of the form 

= {X{r,s,t)/r = 6K{s,t), s G U} (31) 

which is (formally) expected to be a representation of the level set at time 

t of the solution of Problem (P*^). 

The stretched variable p is then defined by 

p = p^(x, t) = di^^t)-^h.{S{x,t),t) ^32) 

which represents the distance from x to 7^ in the normal direction divided 
by e. From now on, we use {p, s, t) as independent variables for the inner 
expansions. The relation between the old and new variables are 

X = X{p,s,t) = X{e{p + he{s,t)),s,t) 

= Xo{s,t)+e{p + h,{s,t)) n{s,t) (33) 

We associate to any function w{x,t) the function 

w{p, s, t) = w{Xo{s, t) + e{p + K{s, t))n(s, t), t) (34) 

or equivalently 

^^ d{x,t) - 6K{S{x,t),t) 
w{x, t) = w{ , S{x, t),t). 
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Note that 

w{r, s, tj = w{ , s, t). 

By differentiation we obtain the following formulas 

dtw = {—Ve~^ — dfhe)wp + d^w 

Vw = (ne^^ - V^he)wp + V^w 

Aw = {e^^ + iV^/ieHwpp + {Ade^^ - A^he)wp 

-2S/^he.V^Wp + A^w, (35) 

where in the above formula for Aw, 

Ad = Ad\x=Xois,t)+e{p+hc{s,t))n{s,t) 

^ K{s,t) - eip + h,is,t))b{s,t) + J2^Ms,t)ip + his,t))\ (36) 

i>2 

with b defined in (l26l) . K defined in ( [25|) and for some given functions 
{bi{s,t))i>2 only depending on 7^. Therefore 

e^{dtw — Aw) = —Wpp — e{V + Ad)wp 

+ e^d'^w - A^w) - (9[/i, - A^h,)wp] 

- e^[\V^K\^Wpp-2V^K.V^Wp] (37) 

The Jacobi For later purposes, we need to compute the Jacobi of the 
transformation X. In the (p, s) coordinates, dx = eJ^{p, s,t)dsdp where ds 
is the surface element of •jt and where eJ'^{p, s, t) = dX{p, s, t)/d{p, s) is the 
Jacobi. We prove below that 

Lemma 4 For all p eR, s e U and t G [0, T], 

n-l 

J'ip,s,t)= + e{p + h%s,t))Ki{s,t)]. (38) 

i=l 

Proof. The equality (l38l) is obtained in two steps. First we consider the 
function X = X{r,s,t) defined in (fTSi) . denote its Jacobi by J = J{r,s,t) 
and prove that for all p e K., s G f/ and t G [0, T], 

J'ip, s, t) = J{e{p + K{s, t)), s, t). (39) 
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Second we compute J and show that for all s & U, for all t G [0, T], 

n-1 

J(r,s,t) = nil + rKi{s,t)] 
1=1 

n— 1 

= 1 + Ad{Xois,t),t)r + J2r'j,{s,t), (40) 

i=2 

for some given functions ji depending on 74. Consequently fl38|) follows di- 
rectly from fl5^ and fHO]) . 

In order to establish (139|) . note that by definition (!33|) . 

X(p,s,t) = X(e(p + /i,(s,t)),s,t) 

so that 

dX _ dX 
dp dr 

and for i = 1, ...,n — 1, 

dX _ dX dh' dX 

dsi dsi dsi dr 

Thus for all p G ffi., s G f/ and t G [0, T], 

, ,dX dX dh'dX dX dh' dX^ 
eJ (p, t) =edet [— , — + e— — , ^ + e^^] 

dX dX dX 
which is (1391) . 

In order to establish fHU]) . we consider the Hessian matrix of d on 7^ and 
denote for s G ?7 and t G [0, T] 

A = A(s,t)=D2^(Xo(s,t),t) 

SO that (1241) reads 

A.n(s,t) = 0. (41) 

Moreover, differentiating the identity (!22l) at r = with respect to Sj for 
z = 1, n — 1 yields 

14 



From 

X{r,s,t) = Xo{s,t) + rn{s,t), 
it follows that using (14T!) 

dX 

— = n{s, t) = {In + rA{s, t))(n(s, t)) 
and that, using ( l42l) for i = 1, n — 1, 

^ = ^ + r— = {In + rA{s, t))(^). 
9si dsi dsi ' 

Therefore for all s G t/ and t G [0, T], 

J{r,s,t) = det K-,^,---,^ 

= det [(/„ + rA)(n), (/„ + rA)(-^), (/„ + rA)(— ^)] 

r T . / M -, r dXn dXf) . 

= det J„ + rA{s, t)] det n, — ^, — ^ 

which in view of (fTTl) proves that 

J{r,s,t) = det [/„ + rA(s,t)] 
which yields (140|) . since the eigenvalues of A{s, t) are ki, ■ ■ ■ , 0. 

5 The approximate solution 
5.1 Asymptotic Expansions 

Let k > max{2,n/2} be a fixed integer. In the sequel, we use the sign ^ 
to represent an asymptotic expansion; namely, 0^ !^ ^ means that for 

j>0 

j 

every integer j G 151, 0^ = ^e*^, + 0{l)e^^^ where 0(1) is bounded inde- 

1=0 

pendently of e G (0,1). For example, since / is smooth, for any bounded 
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sequence {6, oq, ai, a2, ...}, we have the asymptotic expansion 

i>0 j>0 i>0 

^ /(6)+£f(6)E^'«. + ^'E^V.(fe,«o,-,a^) (43) 

i>0 i>0 

where for any fixed b, fi{b,ao, ...,ai) is a polynomial in (ao,...,aj) of degree 
< i- 

Outer expansion We expand A'^(t) and u'^{x,t) for |d(x,t)| > 35 by 

X'{t) ^ Ao(t) + e Ai(t) + A2(t) + ■ ■ ■ (44) 
u'{x, t) ^ uf{t) := ±1 + £{<(t) + euf{t) + ■■■}. (45) 
Substituting (jHD and (|15]) into ([6]) gives 

/(^±(t))=eA^(t)+52(t.±)'(t) 

which yields for alH > 0, 

ufit) = {A, - /._i(±l, 4, ■ ■ ■ , uf-i) - «f_2,t}/r(±l) (46) 

where /_i = ■u^'ia = 0, u% = ±1, and /j {i > 0) is as in (143|) . Hence, are 
determined by {Aq, ■ • • , Aj}. 

Inner expansion We shall assume that /i^ has the asymptotic expansion 

eh%s,t) ^6hi{s,t)+e^h2{s,t) + ---, {s,t) e U x [0,T] (47) 

Near the interface, we assume that the function associated to u"^ by fl34l) 
has the asymptotic expansion 

u%p, s, t) ^ 6'o(p) + e {uo{p, s,t)+e ui{p, s,t) -\ }. (48) 

In the sequel, the zero-th order expansion refers to 

{d(x, t), Ao(t), Uo(p, s, t),u^{t)} 

and the i-th order expansion refers to 

{hi{s, t), Xi{t),Ui{p, s, t),uf{t)}. 

We shall use (■ ■ ■)i_i to denote a generic function of (p, s, t) depending only 
on expansions of order < i — 1. 
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Matching condition We impose that for all i G 

V(s,t) eUx [0,T], lim Ui{p,s,t) = uf{t) (49) 

p— ►±oo 

Translation We also impose for all i G N, 

V(s, t) G [/ X [0, T], Mi(0, s, t) = (50) 
to be consistent with the assumption that p = is the 0- level set of u^. 

5.2 The u-equation in the new variables 

The equation ([6]) is 

-f{u) = -e\ut-/\u)-eK{t)- 

In the new variables {p,s,t), using (157|) . it becomes the following equation 
for the function u = associated to by fl34l) . 

-f{u) = Upp + e[{V{s,t) + Ad)up-X,] 

+ ^[\V^K\\pp-2V^K.V^Upl (51) 
where V{s,t) is given by (B7|) and Ac? is expanded from (15^ as 

Ac? ^ (s, t) - E [K^, t) + '5^-1 (P, 5, t)] , (52) 

with depending only on expansions of order < z — 1 (in particular, 
6o{p,s,t) = pb{s,t)). Note that 6i-i{p,s,t) is a polynomial in p of degree 
< i, whose coefficients are polynomial in {hi, hi_i) with (s, t)-dependent 
coefficients. 

5.3 The recursive i-th equations 

The zeroth order expansion Since do defined in f|T2|) satisfies 
-f{0o) = {Oo)pp, eo(±oo) = ±l, ^^o(0)=0, 

the equation flSTl) is satisfied at zeroth order as well as the matching and 
translation condition fH^ - flSU]) . 
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The first order expansion At first order (e^), tlie equation (IHTj) imposes 
Cuo = {K{s,t) + V{s,t)){eonp) - Ao(t), (53) 
witli C defined in f|T5|) . The solvabihty condition stated in Lemma [3] reads 
{K{s, t) + V{s, t)) f %f{z)dz = 2Ao(t) 



which reads in view of ( 13. ip 

V{s, t) = -K{s, t) + aAo(t) for s e [/ (54) 

which implies in view of ( 127|1 that 

rff = A(i — o"Ao(t) on 7^. (55) 

Moreover equation (J53il has then a unique solution satisfying (H9ll - (j50l) which 
is given by 

Mo(p,s,t) = -Ao(t)^i(p) (56) 
for all (s, t) eU X [0,T]. Note that for all non-negative 171,11,1, 

D;^D:Dl[uo{p,s,t) - 4{t)] = Oie-'^^p^) as p ^ ±00. 

Higher order expansion Plugging the expansions (H3l) . (1471) . ( l48l) into 
(!5T|) and using (l54l) and (!52l) leads to the following identity 

-/(^o) - sf'ieo)iJ2 ^'^^) - E ^V.(^o, wo, w.) 

j>0 i>0 

= + e{Y,e\u,),,) + e[(aAo(t) - + ^^-l)>P " E^'^*](57) 

j>0 j>l j>0 

+ e'Y.^\^^ -dl)u,-e(Ze\^^ -dl)hM + eY.^\u,),) (58) 

j>0 i>l j>0 

+ [£2|Vr/ie|V-2e(E^^Vr^,).VS]- (59) 

i>l 

Define the operator Af^ acting on functions h = h{s, t) by 

Af^h := {dlh - A^h - bh) (60) 

We derive below the {i + l)-th order expansion for i > 1 and obtain the 
following result. 
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Lemma 5 At order e^~^^, with i >1, the equation 137]) imposes 

Cu,= U^{h)e'^ - \{t) + 6i2(V^/ii.V^/i,)^^' + i?,_i(p, s, t), (61) 

with Ri-i only depending on expansions of order < i — 1. Besides -Rj_i(p, s, t) 
is a polynomial in p of degree < i (whose coefficients are polynomial in 
{hi, hi_i, ui, and in their derivatives with respect to (p, s, t) ). 

Proof. First note that using flS^ . the coefficient of order in flFTI) is 

{ui)pp + crAo(t)(Mi_i)p - b{s,t)hi{s,t)9'f^ - \i(t) + (■ • 
(M,)pp-6Mo-A,(t) + (-- ■).-!, (62) 

with (...)j_i depending only on expansions of order < i — 1. Moreover in 
view of fl52|) . it is a polynomial in p of degree < i (whose coefficients are 
polynomial in (/ii, mi, tij_i) and in their derivatives with respect 
to (p, s,t)). 

Next, in view of (1471). the coefficient of order e*"*"^ in fl58l) is 



(A^ - df)u,^, + (9[ - Ar)M; + 

(9[-Ar)Mo + (---).-2- (63) 



To obtain the term of order e*"*"^ in fl59l). note that 



j>l i>2 j=l 



i>3 

SO that 



i>3 i>0 



Hence the coefficient of order e^'^^ in e'^\V^h^\'^Upp is 

fei,2(v^/ii.v%)^;' + (---)^-2 

with 6i^2 = 1 or 2 for i = 1 or i > 2 respectively. 
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Similarly, the coefficient of order e*"*"^ in the term e^'W^h^.V^Up is 

V^/li_l.V^(Mo)p + V^/li_2.V^(Mi)p + ... + V^/li.V^(Mi_2)p 

where the ffist term cancels out since V'"(mo)' = in view of fl56p : thus it only 
depends on expansions of order < i — 2 and appears below in the remainder. 
Finally at order e*"*"^, with i > 1, the equation ( ISTI) reads 

-f'i0o)ui - fi-i{0o,uo, --yUi^i) = {ui)pp-Xi(t) 

with Ri-i only depending on expansions of order < Moreover i?j_i(p, s, t) 
is a polynomial in p of degree < z as described in Lemma [5l 

The solvability condition According to Lemma [3], the equation fl6T]) has 
a solution if and only if the following solvability condition is satisfied. 

V(s, t)eUx[0,T], f C Ui{p, s, t)9'o{p)dp = 0. (64) 

Note that 

/ h2{V^h.V''h,)9M9',{p)dp = bu{V^h.V^h,){s,t) f h{9',f]\p)dp = 
SO that the condition fl64|) reads 

M^{h,) = a\{t)^r,_^{s,t), (65) 

with 

ri_i(s,t) = -- / R^^i{p,s,t){9Q)\p)dp 

I JR 

only depending on expansions of order < — We summarize these results 
in the next lemma. 

Lemma 6 Let k > 1 be given. Assume that for all i < k — 1, ( fgij) has a 
solution Ui satisfying 

n^D';;Dl[ui{p, s, t) - uf{t)] = ©(p^e-^l^l) as p ±oo. (66) 

Also assume that for i = k, {hi{s,t), Xi{t)} satisfies (E^j- Then for i = k, 
(Ely admits a unique solution satisfying Ui{0, s,t) = and (6^) . 
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The proof follows from Lemma [3] and an induction argument and is omitted. 
Just note that in the limit p ±00, the equation = e'^{ul — Au^) + f{u^) — 
^'^^\x=x{p St) becomes the outer expansion equation, so that Mj(±oo,s,t) = 
uf{t). Furthermore since Ri-i is a polynomial in p of degree < fl^B]) is 
satisfied for each i > and (s, t) E U x [0, T]. 



5.4 Equation for A*^. 

To find A^(t), we use an asymptotic expansion for = j^ul{x,t)dx. We 
denote by fi^(t) the two domains separated by '^l defined in (EI]), with = 
dD.j (t) . Hence in view of (!32|) 



fi+(t) = {xEn \ d{x, t) > 35} U {xe Vis I [d{x, t) - eh'{S{x, t),t]> 0} 
= {xen \ d{x, t) > 35} U {xe V^s \ P%x, t) > 0} (67) 

and 



n;it) = n\nt{t) 

= {x en\d{x,t) < -35} U {x eV^sl p%x,t) <0} (68) 

We write 

/ ul{x,t)dx = / ul{x,t)dx + / Uf{x,t)dx (69) 

Jn J\d(x,t)\>3S J\d(x,t)\<35 

where 

/ ul{x,t)dx = / u1{x,t)dx+ / ul{x,t)dx (70) 

J\d{x,t)\<35 ApHx,t)\>^ J|p=(x,t)|<| 

In the sequel we choose < e < Eq small enough so that 

V£e(0,£ol, max \eh%s,t)\ < - (71) 
seu,te[o,T) ' " 2 



Then it follows that 



5 5 
\p'ix,t)\ > - ^ \dix,t)\ > -. 
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Thus if \d{x,t)\ > 35 or \p^{x,t)\ > |, then \d{x^t)\ > | so that at these 
points {x, t), 

ul{x,t) ^ {uiy{t)x{d{x,t)>o} + {u'^y{t)x{d{x,t)<o} 

(exponentially small terms of order 0(e^^) do not affect the asymptotic 
expansion in the e power series). Therefore in view of fl69|) - fl70|) 



l{x,t)dx ^ / [{u%y{t)x{d{x,t)>o} + {u'_y{t)X{d{x,t)<o}] dx (72) 



+ 



|P^(^,i)l<: 



- «)'(i)XMx,t)>o} - {uty{t)x{d{x,t)<o}] dx 



h+ I 



|p^(^,i)l<; 



K - {u\y {t)x{d{x,t)>Q} - {u^-y{t)X{d{x,t)<Q}\dx, (73) 



where 



h = {utymntm + {u:y{m-{t)\. (74) 

In the second integral, we make the change of variables given in fl33|) and 
substitute the expression of in formula (l35l) to obtain 



p|<5/e 



< - {u%y {t)x{d(x,t)>o} - {u'-y {t)x{d(x,t)<o}\ dx 



+ 



0<p<5/£ 

5/e<p<0 

\p\<S/e 



dl[u^{p, s,t) — {u'^_^)(t)]eJ^{p, s,t) dp ds 



dl[u^{p,s,t) - (ui)(t)]eJ'(p, s,t)dpds 
+ (-Ve-^ -dlh,)—eJ'(p,s,t)dpds 

J\o\<S/e OP 



(75) 



Finally, 
where 

l2 = 

and 



p\<5/£ 



/ (x, t)dx ^ Ii + I2 + I3, 
Jn 



df[u'{p, s,t) - «)(t)x{p>o} - (wi)(i)X{p<o}]e^'(p, s,t)dpds 

(76) 

h = / {-Ve-^ -d',K)—er{p,s,t)dpds. (77) 

J\p\<&le dp 
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The calculation for Ii. The boundary of ^j{t) is which according to fl31l) 
is given in local coordinates (r, s) by r = eh'^{s, t). Therefore in view of pUj) . 
we have that 

\n7{t)\ = \Q-{t)\+ f ^'''^ J{r,s,t)drds 

Ju Jo 

« \n-{t)\+J2e'{f h,is,t)ds+{...),^,}, 
where only depends on expansions of order < -j — 1. Hence 

= M-\^-{t)\ 

^ \n+it)\-J2e'{[ h,{s,t)ds + i...),.^}. 
From the outer expansion, 

i>0 i>l 

with {uf_iy{t) given by (H^ and depending only on expansions of order 
<i — 1. Therefore 

where depends only on expansions of order <i — 1. 

The calculation for l2- Using the expression for d^u'^ in formula (I3UI) and 

fl66|) . we compute 

a[[M^(p,s,t) - «)(t)X{p>0} - (w-)(i)X{p<o}] 

i>l 

^eY,^\dt + Stdsj)[ui{p,s,t) - ut{t)x{p>o} - u^{t)x{p<o}] 
i>i j=i 

i>2 

with 0(p*~^e~"'^') depending only on expansions of order < i — 1. Therefore 
by definition of I2 in (!76l) . 

i>3 
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where (...)i-2 depends only on expansions of order < i — 2. 
The calculation for I^. Using the expansions 

dp ''^"h 
{-V - ed'^K) = dt{Xo{s,t),t) -Y^e^dlh 



i>l 



and rewriting the expression of in ( l38l) as 



n-l 



J^(p,s,t) = X{[l + e{p + h'{s,t))K,{s,t)] 

i=l 

^ 1 + At/(Xo(s, t),t)e{p + h^s, t)) + 



i>2 

with depending only on expansions of order < i — 1, we obtain that 

^^v-Edlh,)—r{p,s,t)^ 

dt{Xo{s, t),t)e'o{p) + J2 e'e'o{p){-dfh, + dt{Xo{s, t), t)/i,Ad) + J2 ^'i-h-i 

i>l i>l 

SO that 

h ~ J^J^{e',dt{s,t) + J:,y,e%i-dlh, + dt{s,t)Ad{s,t)h,) + {...),_,]}dpds 
^ 2 / dt{s,t)ds + Y.e'{2 I {-dlK + {d'Ad)hi]ds + {...),^i]. 

Finally, substituting d^ and dfhi by (!55|) and (!65l) . and using /SFhids = 0, 
we obtain 

I [ ut^ f [Ad - aXo) ds + Y.e'{ f [{-b + dt Ad) hi - a X,]ds + {...), ^i} 



Thus the condition ufdx ^ is equivalent to 



aXoit) = Ad(-,t), (78) 
aXi{t) = -[b{-,t) - dt{-,t)Ad{-,t)M-,t) + A,_i(t), z>l (79) 
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where Ai_i{t) depends only on expansions of order < i — 1, and : = 
1^ J^ (p, the average of over 7t parametrized by U. 
Hence, we obtain closed systems for d,hi,- ■ hi, namely 

dt{s,t) = Ad(s,t) - Ad(s,t), (80) 
dfhi = A^hi + bhi - [b{-,t) - dt(-,t)A(-,t)]/ii(-,t) + Ai_i(t) (81) 

onU X [0,T]. 

5.5 Construction of Expansions of Each Order 

We can now use induction to construct each order of expansion as follows: 

1) Zeroth order. Given a smooth initial interface 70, (IHO!) is equivalent to 
the volume preserving mean curvature flow ([5]). By the result established in 
[13J, there is a time T > such that there is a unique smooth solution on 
a time interval [0,T]. Consequently, T = Uo<t<r(7t ^ {t}) the modified 
distance function d are well defined. Set Xo{t) as in (1751) . Mo(p, s, t) as in flS^ 
and u^{t) = Ao//'(±l) as in (HGl) . We obtain the zeroth order expansion 
{d(x, t), Ao(t), Mo(p, s, t), M^(t)}. 

2) Higher order expansion. Fix z > 1. Assume that all expansions of 
order < i — 1 are constructed. Then Aj_i(-) in flSTl) is known. Since jt is a 
smooth hypersurface without boundary, it follows from standard parabolic 
PDE theory [18] that (ISTj) admits a unique smooth solution (assuming an 
initial condition such as hi{-,0) = on ?7 is given). Consequently, we can 
define Xi{t) as in fl79l) . uf as in fH6|) and as the solution of fl6T]) given by 
LemmaO This gives the z-th order expansion {hi{s, t), Aj(t), Uj(p, s, t),uf{t)} 
and completes the induction. 

5.6 Construction of the Approximate Solution 

Now fix an arbitrary positive integer k. We construct an approximate solu- 
tion ti| such that Lemma [2] can be applied. 

Let (5 > be a small fixed constant such that (i) d(x, t) is smooth in a 35- 
neighborhood of F, and (ii) for each t E [0, T], 7^ is at least 36 distance away 



25 



from dQ. We define 

plix,t) = £-Hd(x,t)-Sf+V/i,(S(x,t),t)}, 

fc+1 

4%^,^) = ^o(Pfc) +^E^'«i(Pfc(2^'^)'S(a;,t),t), 

We note that pl,ul\ are smooth in a 35 neighborhood of T. 

Now let ({s) G C°°(K) be a cut-off function (depending only on 6) satisfying 

C(s) = 1 if |s| < 6, C{s) = if |s| > 25, 
< sC'(s) < 4 if (5 < |s| < 26. 

We define the needed approximation solution n| by 

ul{x,t) := C(dK"+[l-C(d)]{<t,+X{d>o}+<t,~X{d<o}}' 
ulix.t) := n^,(x,t) + / K(.,0)-n^,(.,t)} 

for all {x,t) G X [0,T]. Then by construction u| is an approximation of 
order k as needed in Lemma [2j Here we just remark that (i) in the set 
{{x,t) I 6 < ±d{x,t) < 26}, the limiting behavior flBBl) guarantees that 
u%{x,t) = Me^fc* + (t) + 0(e~'^^/'^^^)), vahd also for differentiation, (ii) dnU% = 
on dVtrp since is a function of t near (9^2^; and (iii) the correction 



{M|(.,0)-M|(.,t)} = - / / {ulUy,T)dTdy 

Jn J[o,t] 

is of order 0{e^~^^), valid also for differentiation. 

This completes the construction of the approximating solution, and also the 
proof of Theorem [H 



Appendix A : Proof of Lemma [H 

We first consider the case n > 4 so that p = A/n. The Gagliardo-Nirenberg- 
Sobolev inequality (see [8], Theorem 2, p. 265) states that there exists C > 
such that for every R G H^{fl), 
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2n , , _ 

with 2* = -. Using Poincare-Wirtinger inequality (see |8j, Theorem 1, 

p. 275), it follows that there exists C > such that for every R G H^{Q) with 
J^Rdx = 0, 

\\R\\l^* < C\\VR\\l^. (82) 
Writing Holder inequality, we have that 

\\RfX= I \R?\R\'<{I \R\'^f'^{t \R\^^'f'^' 
Jo, Jn Jn 



and we choose 



to obtain 



^ n-2 2 ' ^ 2 



I <r II R||2 II p||P 



Combined with fl82p . this yields the inequality 

\\RfX < C\\R\\l,\\VR\\h, 
which is the conclusion of Lemma [H 

Next we consider the case that 1 < n < 3 so that p = 1. Schwarz's inequality 
then gives that 

Il^lli3= / \Rm<\\R\\U\R\\L^ 
Jn 

For n = 1,2,3, by Sobolev's imbedding theorem, C L"^, so that there 
exists C > such that for every R G H^{Q), 

\\R\\l^ < C\\R\\hi. 

Using again Poincare-Wirtinger inequality, we finally deduce that there exists 
C > such that for every R G H^{Q) with J^Rdx = 0, 

which concludes the proof of Lemma [H 
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